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ABSTRACT: We present a measurement of rj meson production in photon-photon inter- 
actions produced by electron-positron beams colliding with y/s = 1 GeV. The measure- 
ment is done with the KLOE detector at the (^-factory DA<£NE with an integrated lu- 
minosity of 0.24 ftr 1 . The e+e~ -> e+e" rj cross section is measured without detecting 
the outgoing electron and positron, selecting the decays rj — > tt + tt~tt° and rj — > 7r°7r°7r . 
The most relevant background is due to e + e~ — > 777 when the monochromatic photon 
escapes detection. The cross section for this process is measured as a{e + e~ — > 777) = 
(856 ± 

^stat ^ l^syst) pb. The combined result for the e + e — > e + e rj cross section is 
a{e + e~ — > e + e~rj) = (32.72 ± l-27 s ^ a ^ ± 0.70 S y S ^) pb. From this we derive the partial 
width T(rj — > 77) = (520±20 s ^ a ^ ±13 S y S ^) eV. This is in agreement with the world average 
and is the most precise measurement to date. 

KEYWORDS: e + e~ collisions, rj meson, decay width 
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1 Introduction 

Photon-photon production of neutral mesons provides basic information on their structure. 
The strength of the coupling, measured by the partial decay width T(X — > 77), is related to 
the quark content of the meson and gives information on the relations between the hadronic 
state and its qq representation. For the light pseudoscalar mesons 7r°, 77 and 77', the coupling 
to real photons is measured in their 77 decays, while the coupling to space-like photons 
can be measured in 77 interactions. This is of particular interest in evaluating the light- 
by-light contribution to the anomalous magnetic moment of the muon [1]. Photon-photon 
interactions in electron-positron colliders were pioneered at the Frascati collider Adone in 
the '70s [2-4] and since then have been used to study the production of hadrons in almost all 
e + e~ colliders in a variety of conditions in low- and high-g 2 processes [5-7]. In particular, 
measurements of the 77 partial width of 77 and 77' mesons have been done measuring the 
e + e~ — > e + e~r](r]') cross section [8-13]. 
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We present a measurement of the cross section e + e~ — > e + e~r/ with the KLOE detector 
at the ^-factory DA^NE. The cross section a(e + e~ — > e + e~r]) is a convolution of the 
differential 77 luminosity and the 77—7-7/ cross section. The 77 partial decay width T(r] —> 
77) is obtained by extrapolating the value of 17(77 ~~ ^ r /) f° r rea l photons. 

DA<E>NE is an e + e~ collider designed to operate at high luminosity at the mass of the 
(j> resonance, 1020 MeV. We analyzed data collected with DA<£NE operating off the eft peak, 
at -^/s = 1 GeV, to reduce the large background from (j> decays. The final state e + and 
e~ are not detected, being emitted with high probability in the forward directions outside 
the acceptance of the detector. The production of the i] meson is identified in two decay 
modes, rj — > 7r + 7r~7r and 7/ — > 7r°7r 7r , that exploit in a complementary way the tracking 
system and calorimeter of the detector. The most relevant background is the radiative 
process e + e~ — > 777 and, in both measurements, the yield of 77 mesons is controlled by the 
e + e~ —> 777 cross section measured in the same data sample with a dedicated analysis. The 
data sample used in the analyses corresponds to an integrated luminosity of 0.24 fb _1 . 



2 Signal and background model 

For electron and positron beams colliding with energy E, the cross section for production 
of a state X in 77 interactions with photon 4-momenta q± and q2 is 

CT 77 _>x(<7l,<72) $(gi,<?2) ^T^T> ( 2A ) 

where the 77 differential luminosity ^((71,(72) has been calculated in [14-16] using different 
approximations and is proportional to (a / '2ir) 2 (In E / m e ) 2 . For a narrow resonance of spin 
the formation cross section is 

<7 77 _> x = — r x ^ 77 5(w 2 -m 2 x ) \F(q 2 ,q 2 )\ 2 , (2.2) 
mx 

where rx-5.77 is the radiative width, and w 2 = (q\ + (j^) 2 - The transition form factor, 
F(qf, q%), is equal to one for real photons and is usually parametrized in the form 

inspired by the Vector Dominance Model [17]. The parameter b for the r\ meson has been 
measured at high q 2 values in 77 experiments with single-tagging [18-20] and in the rj 
leptonic radiative decays rj — > [21-23] at low q 2 values, closer to those of this mea- 

surement. The results do not show appreciable dependence on q 2 and the value assumed in 
this analysis, b v = (1.94 ± 0.15) GeV~ 2 , was obtained as an average of the measurements 
at low q 2 . 

The detector response for signal and background events is fully simulated with the 
Monte Carlo (MC) program Geanfi [24]. While Geanfi contains the event generator for 
all background processes, a new generator for e + e~ — > e + e~X events is developed and 
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interfaced to the detector simulation. Events are generated with exact matrix element 
according to full 3-body phase space distributions [25]. This results in the production of 
77 mesons with non negligible transverse momentum. The relative error due to high-order 
radiative corrections to equation (2.1) is estimated to be 1% [26]. All background processes 
have been extensively studied in other analyses. A source of irreducible background is the 
reaction e + e~ — > 777 when the monochromatic photon is emitted at small angles and is 
not detected. The cross section for this process is measured in the same data sample with 
two independent methods and the results agree with each other providing an important 
consistency check of the analysis. The beam-induced backgrounds were measured during 
data taking and background events are added to simulated events in the MC on a run-by-run 
basis. 

3 The KLOE detector 

The KLOE detector consists of a large volume cylindrical drift chamber, surrounded by a 
lead-scintillating fibers finely segmented calorimeter. A superconducting coil around the 
calorimeter provides a 0.52 T axial magnetic field. The beam pipe at the interaction region 
is spherical in shape with 10 cm radius, it is made of a Beryllium-Aluminum alloy of 0.5 
mm thickness. Low-beta quadrupoles are located at ~ 50 cm distance from the interaction 
region. Two small lead-scintillating tiles calorimeters (QCAL) [27] are wrapped around the 
quadrupoles. 

The drift chamber (DC) [28], 4 m in diameter and 3.3 m long, has 12,582 drift cells 
arranged in 58 concentric rings with alternated stereo angles and is filled with a low-density 
gas mixture of 90% Helium-10% isobutane. The chamber shell is made of carbon fiber- 
epoxy composite with an internal wall of 1.1 mm thickness at 25 cm radius. The spatial 
resolutions are a xy ~ 150 /im and a z ~ 2 mm . The momentum resolution for long tracks 
is ct(pt)/pt ~ 0.4%. Vertices are reconstructed with a spatial resolution of ~ 3 mm. 

The calorimeter [29] is divided into a barrel and two end-caps and covers 98% of the 
solid angle. The readout granularity is (4.4 x 4.4) cm 2 , for a total of 2440 cells arranged in 
five layers. Each cell is read out at both ends by photomultipliers. The energy deposits are 
obtained from the signal amplitude while the arrival times and the position along the fibers 
are obtained from the time differences. Cells close in time and space are grouped into energy 
clusters. The cluster energy E is the sum of the cell energies. The cluster time t and position 
f are energy-weighted averages. Energy and time resolutions are (Je/E = 0.057/ \J E (GeV) 
and at = 57 ps/y/E (GeV) © 100 ps, respectively. The cluster space resolution is ctn = 1.4 
cm/ y/E (GeV) along the fibers and <7j_ = 1.3 cm in the orthogonal direction. 

The trigger [30] uses both calorimeter and chamber information. For this analysis the 
events are selected by the calorimeter trigger, requiring two energy deposits with E > 50 
MeV in the barrel or E > 150 MeV in the end-caps. A higher- level cosmic-ray veto rejects 
events with at least two energy deposits above 30 MeV in the outermost calorimeter layer. 

1 KLOE uses a coordinate system where 2 is the bisector of the electron and positron beams, x and y 
define the transverse plane. 
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Data are then analyzed by an event classification filter [24], which selects and streams 
various categories of events in different output files. 

4 Data sample and event preselection 

The data were collected at y/s = 1000.1 MeV with electron and positron beams colliding at 
a small angle with an average transverse momentum of 12.7 MeV in the horizontal plane. 
The average instantaneous luminosity was 7 x 10 31 cm _2 s and the analysis is based on 
an integrated luminosity of 242.5 pb _1 measured with a precision of 0.3% recording large 
angle Bhabha scattering events [31]. 

Data are selected with a background rejection filter [24] before event reconstruction. 
A 1/20 sample of unfiltered data, corresponding to about 11 pb -1 , is also reconstructed 
to define the preselection filter used for the analysis and to evaluate its efficiency for event 
selection. The preselection filter requires 

• at least two energy clusters, neutral (not associated to any track) and prompt (with 
\t - r/c\ < hat); 

• all prompt neutral clusters are required to have energy > 15 MeV and polar angle 
20° < y < 160°; 

• at least one prompt neutral cluster with energy greater than 50 MeV; 

• a ratio of the two highest energy neutral prompt clusters to the total calorimeter 
energy R = (E^ + E l2 )/E tot > 0.3; 

• 100 MeV < Etot < 900 MeV, to reject low energy background events and the high 
rate processes e + e~ — > e + e~(7), e + e~ — > 77. 

5 Cross section for e + e~ — > e + e~rj with r\ — > tt + it~it 
5.1 Event selection 

In addition to the preselection, candidate decays rj — > 7r + 7r~7r should fulfill the following 
requirements 

• two and only two neutral prompt clusters with \t — r/c\ < 3at and polar angle 23° < 
7 < 157°; 

• at least two tracks with opposite curvature that are extrapolated inside a cylinder 
p = \J x 1 + y 2 < 8 cm and \z\ < 8 cm centered around the average beam collision 
point; 

• the distance of the first DC hit to the average beam collision point to be less than 50 
cm for both tracks (in case of two or more tracks with the same curvature, the track 
with best quality parameters is chosen); 

• sum of the two tracks momenta \pi \ + \p2\ < 700 MeV. 
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Figure 1. Correlation between the rj longitudinal momentum, pLrp and the squared missing mass, 
TO mis f° r the reconstructed events that pass the selection cuts of the analysis, for MC signal events 
(left) and data (right). The e + e~ — > 777 events, when the monochromatic photon escapes detection, 
are clearly visible in the data. 



To minimize any selection bias and to optimize the selection efficiency, there is no require- 
ment for the tracks to be associated to clusters in the calorimeter nor that they form a 
vertex. The number of selected events is 3.9 x 10 6 . A small fraction of fully neutral final 
states can survive the two tracks requirement in case of photon conversion 7./V — > e + e~N 
or 7T° Dalitz decay. 

Many background contributions have been considered. The e + e~ — > 777 process is 
a source of irreducible background when r\ decays to 7r + 7r~7r and the monochromatic 
photon, E-f = 350 MeV, is emitted at small polar angles and is not detected. However, the 
correlation of the squared missing mass, m^ is and the 77 longitudinal momentum can be 
used to separate the signal from the background. For the background pl v = E~ cos# ~ 350 
MeV and ~ while the signal, for small values of pT-q, 1S characterized by m^ is ~ 

(\/s — m v ) 2 + {yfs/m^) p\ , as shown in Figure 1. 

The process e + e~ — > ujtt , with co — > 7r + 7r~7r°, has four photons in the final state and 
therefore produces the same final state as the signal when two photons are not detected. 
The cross section has been measured with data from the same run [32], a(e + e~ — > um° — > 
7r + 7r~7T 7r°) = (5.72 ± 0.05) nb. The e + e~ — > KlK$ events can mimic the signal either 
when the Kl decays to ir ± £ Zf i> close to the collision point and K$ decays to 7r°7r°, or when 
the Kl escapes detection and K$ —> 7r°7r° is followed by photon conversion in an e + e~ pair 
or by a n° Dalitz decay. The e + e~ — > K + K~ events can mimic the signal when both kaons 
decay close to the collision point, either — > 7r ± 7r° or — > ir®^v in coincidence with 

— > ^v. Also Bhabha radiative events, e + e~ — > e + e~7, given the large cross section, 
can be a source of background in case of accidental or split clusters. 
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5.2 Reconstruction of r\ — > tt + it tt° decay 

To identify the 7T° meson, clusters are paired choosing the combination that minimizes the 
difference between the two-cluster invariant mass and the tt° mass. This is performed using 
a pseudo-x 2 variable 

2 _ (m 77 - 77V)) 2 .., myy ftTBhyi (TEyj\ ( , 

X ^ ~ ^ Wlth 2 V S~ + ~E~) ' { } 

The energy resolution function is given in Section 3, the 77 invariant mass resolution 
is dominated by the calorimeter energy resolution while the angle measurement gives a 
negligible contribution. Figure 2 shows the distribution of the % 2 7 variable for MC signal 
events and for data. In the following analysis we select events with x 77 < 8. 



1 10< !l 

5 10 5 r 



10 4 r 



5 10 IS 20 25 30 

/'(data) 




X 2 yv (MC signal) 



Figure 2. Distribution of for data (top) and MC signal events (bottom). 

The two tracks momenta are combined with the 7r° to identify r] — > 7r + 7r~7r° decay 
candidates, assigning the charged pion mass to the tracks. A kinematic fit is done requiring 
the invariant mass of 7^77 equal to the 77 mass. In the fit the energies, Ei, the times, tj, 
and the coordinates of the cluster centroid position Xi,yi,Zi, for the two clusters are varied. 
The track momenta are not varied in the minimization since they are measured with much 
better precision than the cluster energies. There are four constraints: the promptness of the 
two clusters, — r^/c = 0, and the mass values m 77 = m^o, m n7T ^ = m„. Figure 3 shows 
the distribution of x 2 from the kinematic fit for MC signal events and for data. We require 
X 2 < 20 to reduce the rj(— > 7r + 7r _ 7r°)7 background. This process has a long tail in the x 2 
distribution due to events with the monochromatic photon in the detector acceptance and 
one photon from the 7T° decay undetected, that are not rejected by the x 2 7 < 8 requirement. 

At this stage of the selection, radiative Bhabha scattering, e + e~ — > e + e~"f, and e + e~ —> 
77 annihilation followed by photon conversion are still a source of background. Separation 
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Figure 3. Distribution of the x% of the kinematic fit for data (top) and MC signal events (bottom). 

of charged pion from electron/positron tracks is done using a likelihood method when a 
cluster is associated to the track [33]. A cluster is associated if the distance between the 
centroid and the extrapolation of the track to the calorimeter wall is less than 50 cm. The 
e-ir likelihood is based on three variables: i) the difference of time of flight, ii) the energy 
of the cluster; Hi) the fractions of energy deposited in the first and in the fifth calorimeter 
layers. In this analysis events with a cluster associated to each track and a value of the 
likelihood estimator log C^/Ce < for both clusters are rejected. The background from 
e + e~ — > 77 annihilation is reduced requiring that the most energetic cluster satisfies the 
conditions E yl < 230 MeV and 27.5° < 6> 7 i < 152.5°. 

Opposite curvature track pairs can originate from split tracks. This is due to the track 
finding algorithm that looks for secondary vertices of kaon decays. Background of split 
tracks is reduced applying a topological cut based on the correlation between the tracks 
opening angle, a nn , and the distance between the first DC hits associated to the two tracks 
by the tracking algorithm. The background from kaon decays in e + e _ — > K^Kg and 
e + e~ — > K + K~ is reduced applying a cut on the two tracks opening angle a nn > 50°. 

Kaon decays are characterized by non prompt energy clusters, thus both the time and 
the energy assigned by the kinematical fit to the neutral clusters are modified by the fit 
constraints. This effect is observed in the time and energy pulls built with the two neutral 
clusters 

,2 _ \^ 



j.meas\2 I jrfit Tpmeas\2 



at — Oi 

i€27 * i€27 i 



where the superscript meas and fit indicate the values measured and returned by the fit, 
respectively. The pulls are required to satisfy < 7 and x% < 8. 
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Final state 


Selection 


vjiooai 




efficiency (7oj 


emciency (7oJ 


e T e rj 


odd i o 

34.4 ± 0.3 


20.8 ± 0.3 


Tl(—> TT^TT 7T U J7 


13.3 


1.93 


7j(— > 7T + 7T~7)7 


/ion 
43.9 


U.UU9U 


^(—t- neutral)7 


0.185 


0.00030 




3.08 


0.023 


KlKs 


0.169 


0.0059 


K+K- 


0.423 


0.0075 


e + e~7 


0.447 


< 0.0004 



Table 1. Selection efficiency, in %, for the signal and the most relevant backgrounds. The column 
Selection includes the efficiency of the trigger, the background filter and the data filters described in 
Sections 4 and 5.1. The global efficiency for e + e~7 is derived before the cut on the e-n likelihood. 

The selection efficiencies are evaluated with the MC simulation and are listed in Table 1 
for the signal and the most relevant background sources. The column Selection includes the 
efficiency of the trigger, the background filter and the data filters described in Sections 4 
and 5.1. The trigger efficiency is controlled by comparison of the calorimeter trigger with a 
complementary trigger based on the drift chamber hit patterns [30]. A sample of unfiltered 
data is used to control the filter efficiency. 

The signal is simulated with different values of the b v parameter 2 of the form factor in 
equation (2.3) and the fit to derive the signal yield is repeated for each value. The values 
of efficiencies shown in Table 1 correspond to = 1.94 GeV -2 . 

5.3 Cross section evaluation 

The analysis cuts described in Section 5.2 select 2977 events. The number of signal events is 
derived with a 2-dimensional fit to the data. The variables used to discriminate the signal 
from background are the squared missing mass and the r] transverse momentum in the 
interval -0.15 GeV 2 < m 2 mis < 0.25 GeV 2 and p Tr] < 300 MeV that contains 2720 events. 
The fit to the data is done using the simulated shapes for the signal and backgrounds, the 
weights are left free except for rj(— > 7T + 7r~7r°)7 whose cross section and error, measured in 
the same data sample (see Section 8), is used as a constraint in the fit. The fit returns the 
fraction of data events fi = Ui/ntot with the constraint = 1. 

The projections of the iT? mis x py„ distribution are shown in Figure 4 for the data 
and the backgrounds weighted by their fractions fi, and the distribution is shown in 
Figure 5. The most relevant background is e + e~ — > 777 characterized by ~ and 

PLrj — ±350 MeV. Table 2 lists the fraction of events returned by the fit using the signal 
efficiency evaluated with b v = 1.94 GeV~ 2 , the fit is repeated for all the other b v values used 
in evaluating the efficiency. The distributions of the variables used in the event selection 
are compared for data and MC simulation, weighted by the fractions fi returned by the fit, 
and good agreement is observed. The fit finds 394 ± 29 signal events. 

2 The values chosen are b v = 0, 0.7, 1.0, 1.5, 1.64, 1.80, 1.94, 2.00, 2.24 GeV -2 
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Figure 4. Projections of the 2-dimensional fit. Left: distribution of the transverse momentum 
of the 7r + 7r~77 system. Right: distribution of the squared missing mass. The contribution of the 
signal is blue, e + e~ — > 777 is red, e + e~ — > urn is black, e + e~ — > e + e~7 is green, e + e~ — > K + K~ 
is light blue and e + e~ — > KgKj^ is purple. 




P L (MeV) 



Figure 5. Distribution of the longitudinal 7r + 7r~77 momentum. The contribution of the signal is 
blue, e + e~ — > 777 is red, e + e~ — > um° is black, e + e~ — > e + e~7 is green, e + e~ — > K + K~ is light 
blue and e + e~ — > KsKl is purple. 



The contributions to the systematic error are evaluated by varying the analysis cuts 
by the r.m.s. width of the distributions of each variable: Xyyi x\i Eyii 9yi> a im, Xt> x\i 
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l*inal state 


fraction 01 events (7oJ 


e~*~e — 


14. 4» ± l.UD 


777 


32.02 ± 0.54 




20.48 ± 1.81 


KlKs 


11.36 ± 1.70 


K + K~ 


15.13 ± 1.81 


e + e~7 


7.54 ± 0.87 



Table 2. Fraction of events, in %, for the signal and the most relevant backgrounds. 



Variable 


Range 


5a/a(%) 


Xyy 


6.6 - 10.8 


+0.67 -0.73 


X-q 


18.5 - 23.5 


+0.06 -0.68 


E-yl 


210 MeV - 250 MeV 


-1.17 -0.33 




26. 5 /153. 5° - 28. 5°/151. 5° 


+1.21 +0.46 


Xt 


6 - 8 


+1.10 -1.22 


xl 


7-9 


+1.89 -1.39 


^7+7T — 


48° - 52° 


-0.21 +0.20 



Table 3. Systematic errors determined varying the cuts for each variable for the cr(e + e — > 
e + e~?/ — » e + e _ 7r + 7r _ 7r°) measurement. 

accounting for their correlation. This results in a systematic relative error of —2.4% to 
+2.6%. The contributions are listed in Table 3. 

The MC simulation statistical error of 1.4% (Table 1) is added in quadrature, the 
uncertainties in the form factor and in the branching ratio are kept separate to account for 
correlations between the two 77 decay modes. The change of the result due to the variation 
of b v in the transition form factor formula leads to a 2.0% fractional error. We obtain 
a(e + e~ -»• e + e~r] -+ e + e~ it + it~ tt°) = (7.84 ± 0.57 stat ± 0.23 syst ± 0.16 FF ) pb. Using for 
the branching fraction the value BR(rj -)• 7r + 7r~7r ) = 0.2274 ± 0.0028 [34], we obtain 

a(e + e - ^ e + e - v ) = (34.5 ± 2.5 stat ± 1.0 syst ± 0.7 FF ± 0.4 BR ) pb . (5.3) 

6 Cross section for e + e~ — > e + e~rj with r\ — > 7r 7r°7r° 
6.1 Event selection 

In addition to the preselection described in Section 4, candidate decays rj — > 3tt° should 
fulfill the following requirements 

• six and only six neutral prompt clusters with E~ > 15 MeV, \t — r/c\ < 3<Tt and polar 
angle 23° < y < 157°; 

• no tracks in the drift chamber. 
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Final state 


selection 


Global 




emciency /o j 


emciency (yo J 




e + e 7] 




Z5.D ± 0.3 


T}\—} 7T 7T 7T J7 


to.y 


Z.14 


U>7T° 


0.145 


0.0077 


K L K S 


0.0126 


0.0073 


a (980)7 


2.70 


0.85 


fo (980)7 


0.147 


0.0070 


n'i 


2.13 


0.212 



Table 4. Selection efficiency, in %, for the signal and the most relevant backgrounds. The column 
Selection includes the efficiency of the trigger, the background filter and the data filters described 
in Sections 4 and 6.1. 

The number of selected events is 9857. Many background contributions have been 
considered. As in the charged decay analysis, the e + e~ — > rjj process is a source of irre- 
ducible background when rj decays to 3ir° and the recoil photon is not detected. The process 
e + e~ — > W7T with uj — > 7r°7 produces 5 photons in the final state and is important in case of 
accidental or split clusters. The cross section has been measured with the same data set [32]: 
a{e + e~ -> wvr -> vr 7r°7) = (0.550 ± 0.005) nb. The process e+e" -> a (980)7 -> i]ir°j 
can mimic the signal when r) decays to 3-7T and three photons are not detected, or it de- 
cays to 77 with split or accidental clusters. Similarly for e + e~ — > /o(980)7 — > and 
e + e~ — > T] >r y when rj decays to neutrals. Also the process e + e~ — > K^Kg with Kg — > 7r°7r° 
and undetected Kj_, can mimic the signal in case of split or accidental clusters. 

6.2 Reconstruction of r\ — > 3tt° decay 

The six photons are paired choosing the combination that minimizes the difference between 
the 77 invariant mass of the pairs and the mass of the ir° as described in Section 5.2. In 
the following analysis we select events with Xq^ < 14- A kinematic fit is done requiring the 
67 invariant mass to be equal to the rj mass. In the fit the energies, Ei, the times, tj, and 
the coordinates of the centroid positions Xi,yi,Zi, for the six clusters are varied. There are 
seven constraints: the promptness of the six clusters ti — ri/c = and mg 7 = m„. Figure 6 
shows the distribution of \% from the kinematic fit for MC signal events and for data, we 
require xt < 20 to reduce the background. 

MC simulation shows that e + e~ — > 777 gives a large contribution to the tail of the 
distribution when the monochromatic photon is in the acceptance and is wrongly paired 
with a photon from r\ decay. In this case it also produces an enhancement at large values of 
the 67 invariant mass distribution. To reduce the background we require the highest energy 
neutral cluster to have < 260 MeV and the six-photon invariant mass < 630 MeV. 

The selection efficiencies are evaluated with the MC simulation described in Section 2 
and are listed in Table 4 for the signal and the most relevant background sources. 
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Figure 6. Distribution of for data (top) and MC signal events (bottom). 
6.3 Cross section evaluation 

The number of signal events is derived with a 2-dimensional fit to the data. The distribu- 
tions used to discriminate the signal from background are the squared missing mass and 
the rj longitudinal momentum in the interval -0.15 GeV 2 < m 2 ^ < 0.35 GeV 2 and -450 
MeV < pl v < 450 MeV that contains 2166 events. The fit to the data is done using the 
simulated shapes for the signal and backgrounds and the fit returns the fraction of data 
events /, = rii/ntot with the constraint fa = 1. 

The contribution of all backgrounds, except rj'y production, is very small, below the 
statistical sensitivity of the fit. The contribution of lutt derived from the value of the 
e + e~ — > W7T° cross section [32] is / W7r o = 0.47%. The contributions of ao(980)7, 7/7 and 
KlKs expected extrapolating the measurements at the <f> peak are negligible. 

The fit with two components gives f eer) = (33.4 ± 1.5)% and = (66.6 ± 1.9)% using 
the signal efficiency evaluated with b v = 1.94 GeV -2 . The fit is repeated for all the other 
values. The projections of the m^is x pi v distribution are shown in Figure 7 for the data 
and the background weighted by their relative factors /j, and the pttj distribution is shown 
in Figure 8. 

The contributions to the systematic error are evaluated by varying the analysis cuts 
by the r.m.s. width of the distributions of each variable: X77) Xr)> ^-yli 7Ti 67> accounting 
for their correlation. This results in a systematic relative error of —1.5% to +2.6%. The 
contributions are listed in Table 5. 

The MC simulation statistical error of 1.0% (Table 4) is added in quadrature, the 
errors due to knowledge of the form factor and to the branching ratio are kept separate. 
The changes of the result due to the variation of b v in the transition form factor formula lead 
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Figure 7. Projections of the 2-dimensional fit. Left: distribution of the 67 longitudinal momentum. 
Right: distribution of the squared missing mass. The contribution of the signal is blue, e + e~ — > 777 
is red. 
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Figure 8. Distribution of the 67 transverse momentum. The contribution of the signal is blue, 
e + e~ —¥ 777 is red 

to a 0.7% fractional error. We obtain a(e + e~ — > e + e~rj — > e + e~37r°) = (10.43 ± 0.48 s t a t ± 
0.29 S y S ^ ± 0.07pp) pb. The analysis of the systematic uncertainties of the e + e~ — > 777 
measurement leads to a relative error of 0.6%: we obtain a{e + e~ — > 177 — > 37r°7) = 
(278.0±8.1g^ a ^±1.7 S y S ^) pb. Using for the branching fraction the value BR(r] — > 7r°7r 7r°) = 
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Variable 


Range 


Sa/a(%) 


A-yy 


12 - 16 


-0.51 +0.83 


X.T) 


17 - 23 


-0.26 -0.68 


Af 67 


610 MeV - 650 MeV 


-1.33 +2.38 



Table 5. Systematic errors determined varying the cuts for each variable for the a(e + e — > 
e + e~ri — > e + e~37r°) measurement. Varying the cut on E 7 i gives a negligible contribution. 

0.3257 ± 0.0023 [34], we obtain 

a(e + e~ -> e + e~ V ) = (32.0 ± 1.5 stat ± 0.9 syst ± 0.2 FF ± 0.2 BR ) pb (6.1) 

and 

a{e + e - ^ 7?7 ) = (853 ± 25 stat ± 5 syst ± 6 BR ) pb . (6.2) 
7 Determination of T(rj — > 77) 

The two values of the cross section in equations (5.3) and (6.1) are combined accounting 
for the following sources of correlation: 

• systematic uncertainties are correlated due to the requirements on the neutral prompt 
clusters, the photon energy, time and position resolutions common to both selections 
and fit procedures; 

• the determination of the signal efficiencies for the two measurements that share the 
same transition form factor; 

• the systematic error in the measurement of the luminosity [31]; 

• the correlation between the r\ — > 7r + 7r _ 7r° and 77 — > 3tt° branching ratios [34]. 
From the combination of the two measurements we derive 

a ( e + e - e+e-r,) = (32.7 ± 1.3 stat ± 0.7 syst ) pb . (7.1) 

The partial width of the 77 meson, T(r] — > 77), can be determined from equations (2.1) 
and (2.2). The 77 differential luminosity is calculated following reference [25], the program 
computes also the transition form factor as parametrized in equation (2.3), for the same 
values of the b„ parameter used in evaluating the e + e~ — > e + e~r] cross section. Since the 
values of the 4-momenta qi and 52 sampled in the two decay modes analyzed in Sections 5 
and 6 can be slightly different, the partial width is determined separately for the two decays. 
The theoretical error in evaluating 17(77 — > rj) has been added to the systematic error due 
to the form factor. From the two values of the e + e~ — > e + e~ry cross section, (5.3) and (6.1), 
we derive 

V -> vr+vr-^ T( V -> 77) = (548 ± 40 stat ± 16 syst ± 14 FF ± 7 BR ) eV , 
v ^ vrO^o^o r{rj ^ 77) = (5Q9 ± 23gtat ± 14gyst ± 8pp ± 4br) eV 

The two measurements are combined accounting for their correlations to derive 

rfa -> 77 ) = (520 ± 20 stat ± 13 syst ) eV . (7.3) 
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8 Measurement of the cross section for e + e — > 777 

The most relevant background in the measurement of the e + e~ — > e + e~r] cross section is 
due to the radiative process e + e~ — > 777. The value of the cross section has been used as 
a constraint in the fit in case of the r] — > tt + tt~tt° decay while it has been derived as a 
by-product of the analysis of the 77 — > 3ir° decay. The cross section has been measured by 
the SND experiment [35] at VEPP-2M in the range y/s = (0.6 - 1.38) GeV, but with less 
precision than needed to control the analysis of e + e~ — > e + e~rj. 

The cross section for e + e~ — > 777 is measured exploiting the 77 — > 7r + 7r _ 7r° decay 
using the same data sample and the same preselection procedure described in Sections 4 
and 5.1 with the only difference that in this case events with three and only three neutral 
prompt clusters are selected. The event selection aims at finding two tracks of opposite 
curvature, compatible with being due to n^, two neutral prompt clusters compatible with 
being originated by a 7r° decay, and a third neutral prompt cluster compatible with the 
photon recoiling against the 7r + 7r~7r° system. 

Several background processes have been considered. e + e~ — > ujtt with oj — > 7r + 7r~7r° 
is characterized by two tracks and four photons and can simulate the signal if one photon is 
not detected. e + e~ — > 7r + 7r _ 7r°7 has the same configuration as the signal. e + e~ — > K^Ks 
can mimic the signal when Kl decays to ir^t^v close to the collision point and K$ decays 
to 7T 7r° but one photon is not detected. e + e~ — > K + K~ can mimic the signal when both 
kaons decay close to the collision point to tt^tt , 7r T 7r° and one photon is not detected, or 
decay to tt^tt , and the additional photon originates from split or accidental clusters. 
e + e~ — > tt + tt^tt° and e + e~ — > 7r + 7r~7 can mimic the signal in case of one or two accidental 
or split clusters. e + e~ —> e + e~7 has a very large cross section and can be an important 
background if the electron (positron) is misidentified as a pion and the two additional 
photons originate from split or accidental clusters. e + e~ —> 77 has also a large cross 
section and may originate background in case of photon conversions and there are split or 
accidental clusters. Beside these, 777 production with 77 decaying to 7r + 7r~7 or to 3ir° should 
be discriminated from the 7r + 7r~7r° signal by the number of prompt neutral clusters. 

8.1 Reconstruction of 7/7 — > 7r + 7r~7r°7 events 

The identification of the 7r° meson follows the procedure described in Section 5.2. No cut 
is applied to the value of ^77- A kinematic fit is applied to the selected combination of 
three neutral prompt clusters and two tracks, with the assignment of the charged pion 
mass. The fit uses 15 variables, the energy Ei, time t{ ctnd cluster coordinates xi^yi^Zi of 
the three clusters, and has 7 constraints, promptness of three clusters ij — rj/c = 0, energy 
and momentum conservation: E^i + E n + + E n - = y/s and J2iP"/i +P71— = Pe+e~ ■ 

The track momenta are not varied in the minimization procedure. Figure 9 shows the 
distribution of the x 2 °f the kinematic fit for MC signal events and for data. In the 
following analysis we select events with x 2 < 50. 

The background of e + e~ — > e + e~7 and e + e~ —> 77 is reduced using the e-ir likelihood 
estimator as described in Section 5.2, and requiring the angle between the two tracks to 
be a^,,- < 160°, and the angle between any photon pair to be a 77 > 20°. The background 
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Figure 9. Distribution of x 2 for the data (top) and MC signal events (bottom). 

of e + e~ — )• 77-7 — t- 37r°7 with photon conversion is reduced requiring the sum of the photon 
energies E™ < 660 MeV. At this stage of the analysis, the residual background is 
dominated by the processes e + e~ — > 7r + 7r~7 and e + e~ — > 7r + 7r~7r with split clusters, 
characterized by a neutral energy smaller than for the signal, and e + e~ — > a; 7 — > 7r + 7r~7r°7 
characterized by the same final state as the signal. These backgrounds are reduced by 
requiring for the sum of the track momenta \p + \ + < 440 MeV. The effect of these cuts 
is controlled by the distribution of the energy of the unpaired photon shown in Figure 10 
where is the value returned by the fit and has a resolution greatly improved by the good 
time and position resolution of the calorimeter. The peaks at the energies of the photon 
recoiling against the uj and the rj are clearly visible over a small background at £^3 = 194 
MeV and E^ = 350 MeV, respectively. 

The selection efficiencies are evaluated with the MC simulation described in Section 2 
and are listed in Table 6 for the signal and the most relevant background sources. 

8.2 Evaluation of the cross section 

The number of signal events is derived with a 2-dimensional fit to the data. The dis- 
tributions used to discriminate the signal from background are the energy of the un- 
paired photon and the invariant mass of the two charged pions in the interval 50 MeV 
< £ 7 3 < 400 MeV and 280 MeV < m n7r < 520 MeV that contains 55150 events. The fit 
to the data is done using the simulated shapes for the signal and backgrounds and the 
weights are left free. The projections of the E^ x m n7T distribution are shown in Fig- 
ure 11 for the data and the backgrounds weighted by their relative factors returned by 
the fit. The result of the fit gives 13536 ± 121 signal events resulting in a cross section 
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Figure 10. Distribution of the energy of the unpaired photon for data before (black) and after 
(red) the cut on the sum of the tracks momenta. The e + e~ — > ujj peak is clearly visible, with 
£L = 194 MeV. 



Final state 


Selection 


Global 




efficiency (%) 


efficiency (%) 


r)(— > 7r + 7r~7r°)7 


36.71 ± 0.02 


28.68 ± 0.02 


7r + 7r~7r°7 


6.08 


1.19 


U>7T° 


19.80 


1.07 


T]{— > 7T + 7T~7)7 


0.723 


0.069 


r)(— > neutral)7 


0.111 


0.002 



Table 6. Selection efficiency, in %, for the signal and the most relevant backgrounds. The column 
Selection includes the efficiency of the trigger, the background filter and the data filter described 
in Section 4. 

a(e + e~ -)• 777 vr + 7r-7r 7) = (194.7 ± 1.8 stat ) pb. 

The only relevant backgrounds are from e + e~ — > 7r + 7r _ 7r°7 and e + e~ — > urn — > 
7T + 7r~7r 7r°. The distributions of the signal and e + e~ — > 7r + 7r _ 7r°7 are well reproduced 
both in shape and relative normalization, while the fraction of a>7r° events results slightly 
higher than expected. If the measured value and its error, cr(e + e~ — > cjtt — > 7r + 7r~7r°7r°) = 
(5.72 ± 0.05) nb [32], are introduced as a constraint, the fit returns a value 1.36% higher 
for the e + e~ — > 7/7 cross section. This difference is accounted for in the systematic error. 

Other contributions to the systematic error are evaluated by varying the analysis cuts by 
the r.m.s. width of the distributions of each variable, x 2 , a 77 , + \p~\, accounting 
for their correlation. This results in a relative error of ±1.45% and a{e + e~ — > 777 — > 
7r + 7r~7r°7) = (194.7 ± l-8 s ^ a ^ ± 2.8 S y S ^) pb. Using the branching fraction for r\ — > 7r + 7r _ 7r , 
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Figure 11. Projections of the 2-dimensional fit. Left: distribution of the energy of the unpaired 
photon. Right: distribution of the invariant mass m w + n - . The contribution of the signal e + e~ — > 7/7 
is blue, e + e~ — > lutt is green and e + e~ — > ir + ir~ir°~f is purple. 

we derive 

ff ( e + e - ^ m) = ( 856 ± 8gtat ± i2 sygt ± H BR ) pb (8.1) 

This value, obtained from a direct measurement, agrees well with the value (6.2) obtained 
from the analysis of 77 — > r\ — > 3tt°. The result interpolates well with the measurements of 
the SND experiment [35] and has a better precision. 

9 Summary 

The cross section a(e + e~ — > e + e~r]) has been measured at -^/s = 1 GeV with the KLOE 
detector based on an integrated luminosity of 0.24 fb _1 . The r] mesons are selected using 
the two decays 77 — > 7r + 7r~7r and rj — > ir°ir°iT that exploit in a complementary way the 
tracking and the calorimeter measurements. Many background processes are considered, 
the most relevant being e + e~ — > 777 when the photon is emitted at small polar angles 
and escapes detection. As a consistency check, we have measured the cross section for 
e + e~ — > rj^f in two independent ways, the two values agree well with each other and we 
derive a(e + e~ — > 777) = (856 ± 8 s ^ a ^ ± 16 svs t) pb. This value interpolates well previous 
measurements by the SND experiment and is more precise. The cross section for e + e~ — > 
e + e~rj is obtained independently for the two 77 decay modes with a 2-dimensional fit to the 
squared missing mass and the r/ momentum projections. Combining the two measurements 
we obtain a{e + e~ — > e + e~rf) = (32.72 ± 1.27 s ^ a ^ ± 0.70 S y S ^) pb. This value is used to 
extract the partial width T(r/ — > 77) = (520 ± 20 s ^ a ^ ± 13 S y S ^) eV. This is in agreement 
with the world average of (510 ± 26) eV and is the most precise measurement to date. 
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